Solving Polynomial Systems via a Stabilized Representation of Quotient Algebras

We consider the problem of finding the isolated common roots of a set of polynomial functions defining a zero-dimensional ideal I in a ring R of polynomials over C. We propose a general algebraic framework to find the solutions and to compute the structure of the quotient ring R/I from the null space of a Macaulay-type matrix. The affine dense, affine sparse, homogeneous and multi-homogeneous cases are treated. In the presented framework, the concept of a border basis is generalized by relaxing the conditions on the set of basis elements. This allows for algorithms to adapt the choice of basis in order to enhance the numerical stability. We present such an algorithm and show numerical results

Polynomial eigenvalue solver based on tropically scaled Lagrange linearization

We propose an algorithm to solve polynomial eigenvalue problems via linearization combining several ingredients: a specific choice of linearization, which is constructed using input from tropical algebra and the notion of well-separated tropical roots, an appropriate scaling applied to the linearization and a modified stopping criterion for the $QZ$ iterations that takes advantage of the properties of our scaled linearization. Numerical experiments suggest that our polynomial eigensolver computes all the finite and well-conditioned eigenvalues to high relative accuracy even when they are very different in magnitude.status: publishe

Rank-deficient submatrices of Fourier matrices

AbstractWe consider the maximal rank-deficient submatrices of Fourier matrices with order a power of a prime number. We do this by considering a hierarchical subdivision of these matrices into low rank blocks. We also explore some connections with the fast Fourier transform (FFT), and with an uncertainty principle for Fourier transforms over finite Abelian groups

A generalization of Floater--Hormann interpolants

In this paper the interpolating rational functions introduced by Floater and Hormann are generalized leading to a whole new family of rational functions depending on $\gamma$, an additional positive integer parameter. For $\gamma = 1$, the original Floater--Hormann interpolants are obtained. When $\gamma>1$ we prove that the new rational functions share a lot of the nice properties of the original Floater--Hormann functions. Indeed, for any configuration of nodes, they have no real poles, interpolate the given data, preserve the polynomials up to a certain fixed degree, and have a barycentric-type representation. Moreover, we estimate the associated Lebesgue constants in terms of the minimum ($h^*$) and maximum ($h$) distance between two consecutive nodes. It turns out that, in contrast to the original Floater-Hormann interpolants, for all $\gamma > 1$ we get uniformly bounded Lebesgue constants in the case of equidistant and quasi-equidistant nodes configurations (i.e., when $h\sim h^*$). In such cases, we also estimate the uniform and the pointwise approximation errors for functions having different degree of smoothness. Numerical experiments illustrate the theoretical results and show a better error profile for less smooth functions compared to the original Floater-Hormann interpolants.Comment: 29 page

Suppose I Met You Face To Face

https://digitalcommons.library.umaine.edu/mmb-vp/2936/thumbnail.jp

Bézout and Hankel matrices associated with row reduced matrix polynomials, Barnett-type formulas

AbstractBased on the approach introduced by B.D.O. Anderson and E.I. Jury in 1976, the definition of finite Hankel and Bézout matrices corresponding to matrix polynomials is extended to the case where the denominator of the corresponding rational matrix function is not necessarily monic but is row reduced. The matrices introduced keep most of the well-known properties that hold in the monic case. In particular, we derive extensions of formulas giving a connection with polynomials in the companion matrix (usually called Barnett formulas), of the inversion theorem and of formulas concerning alternating products of Hankel and Bézout matrices

An Approximate Projection onto the Tangent Cone to the Variety of Third-Order Tensors of Bounded Tensor-Train Rank

An approximate projection onto the tangent cone to the variety of third-order tensors of bounded tensor-train rank is proposed and proven to satisfy a better angle condition than the one proposed by Kutschan (2019). Such an approximate projection enables, e.g., to compute gradient-related directions in the tangent cone, as required by algorithms aiming at minimizing a continuously differentiable function on the variety, a problem appearing notably in tensor completion. A numerical experiment is presented which indicates that, in practice, the angle condition satisfied by the proposed approximate projection is better than both the one satisfied by the approximate projection introduced by Kutschan and the proven theoretical bound